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Abstract: Plasmonic dark modes are pure near-field resonances since 

their dipole moments are vanishing in far field. These modes are particularly 
interesting to enhance nonlinear light-matter interaction at the nanometer 
scale because radiative losses are mitigated therefore increasing the intrinsic 
lifetime of the resonances. However, the excitation of dark modes by 
standard far field approaches is generally inefficient because the symmetry 
of the electromagnetic near-field distribution has a poor overlap with the ex¬ 
citation field. Here, we demonstrate the selective optical excitation of bright 
and dark plasmonic modes of single gold nanorods by spatial phase-shaping 
the excitation beam. Using two-photon luminescence measurements, we 
unambiguously identify the symmetry and the order of the emitting modes 
and analyze their angular distribution by Fourier-space imaging. 
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Introduction 

Localized surface plasmons (LSPs) are quantized collective oscillations of conduction electrons 
sustained by metal nanoparticles. The electromagnetic fields associated with LSPs are typically 
enhanced and strongly confined to the surface. Manipulation of the surface plasmons triggered 



a paradigm to control light-matter interactions at the nanoscale by a proper engineering of 
the surface EElElll. The development of optical antennas m is largely relying on such 
engineering with benchmark applications ranging from surface-enhanced spectroscopies m to 
local manipulation Qm, and molecular control EKHUUIia. 

LSP responses are generally inferred from the multipolar nature of the charge oscillations 0] 
[13 . Depending on the symmetry of the modes, LSPs may have a dipole-active character and are 
considered as bright modes because they efficiently decay in the far held. Maxwell’s equations 
also predict the existence of dipolar-inactive LSPs : quadrupolar and higher-order eigenmodes 
cannot interact classically with light and they are generally referred as dark modes ifTSl fT^ . 
An interesting consequence is that the intrinsic lifetime of the mode is typically longer than 
radiative modes, providing thus enhanced near-held interactions with single emitters ifT^ or 
other plasmonic elements. As an example, Solis et al showed that light propagation through 
nanoparticle assemblies via dark LSPs is ten fold more efficient than with bright SPs ca. 

Dark resonances may arise when nanoparticles are paired. Near-held coupling results to an 
hybridization of the system with the appearance of bright bonding modes and a dark antibond¬ 
ing modes in the gap. The spatial distribution of these dark modes were mapped by detecting 
nonlinear photoluminescence with focused helds (161, phase-shaped excitations CtJIISI as well 
as second-harmonic generation oa. At the level of an isolated plasmonic structure, an optical 
excitation of dark LSPs is symmetry-forbidden for a normally incident plane-wave excitation. 
The spectral signature of dark modes was however reported with an oblique incidence (^1^ 
and their spatial distribution interrogated by electron-energy loss spectroscopy and cathodolu- 
minescence experiments (^1^ . 

In this letter, we explore the bright and dark modal distributions sustained by isolated sin¬ 
gle gold nanorods as a function of their length. A two-photon luminescence non linear process 
(TPL) is used to probe the local electromagnetic held enhancement provided by LSP reso¬ 
nances EilEllESlEl. Specihcahy, we demonstrate a selective excitation of bright and dark 
SP resonances using a hrst-order Hermit-Gaussian beam. Experimental data are compared to 
numerically calculated near-held intensities and phase distributions. Finally, we demonstrate 
the dipolar and quadrupolar character of resonant bright and dark LSP modes by measuring the 
angular distribution of the back-scattered light at the excitation wavelength. 

Experimental 

Nanorods were realized by electron beam lithography and gold evaporation on a conductive 
glass cover slip. The dimensions of the nanorods were 35 nm thick, 55 nm wide with length L 
varying from 90 nm to 925 nm by steps of 15 nm. A scanning electron micrograph (SEM) of a 
500 nm long nanorod is displayed in inset of Fig.[2 

TPL confocal images were recorded using the experimental arrangement sketched in Fig. 
The patterned cover slip was placed on a piezo-electric scanning stage htted on an inverted mi¬ 
croscope. A TiiSapphire laser producing 120 fs pulses at a wavelength of 800 nm was focused 
on a diffraction-limited spot by high aperture microscope objective (x60, numerical aperture 
NA=L49). Two excitation conhgurations were used to selectively excite bright and dark LSPs. 
Bright modes were excited by a classical gaussian beam prohle (axis v and y are 

dehned in inset of hgure[^. Dark modes were generated by a HG\^{){x^y) beam prohle formed 
by centering a A/4-step gold mirror in the laser beam prior to the microscope objective (28l. 
The indices of these Hermitt-Gaussian beams correspond to the number of ;r-phase changes 
along the corresponding Cartesian direction. For the remaining of the discussion, the excitation 
polarization direction is hxed along the long nanorod axis {x direction according to the inset of 
Fig.[^. The TPL signal generated when the nanorods are raster scanned through the waist of the 
focused beam was recorded by an amplihed photo-multiplier tube (PMT) via a 200 jim optical 



Fig. 1. Schematic view of the experimental setup described in the text. In inset is shown a 
SEM image of a 500 nm long gold nanorod. 


fiber core. The TPL signal was filtered out with the help of a dichroic beamsplitter. To analyze 
the angular distribution at the excitation wavelength, a conjugate Fourier plane of the micro¬ 
scope Z" was projected on a charge-coupled device (CCD) camera. For these measurements, 
we reduced the excitation numerical aperture (NA^xc < 1.0) by a diaphragm placed in a conju¬ 
gate Fourier plane Z' before the microscope. In the detection path, a solid mask positioned in 
an intermediate Fourier plane Z' rejects the direct refiected light within the to only detect 
the scattering contribution distributed between the NAexc and the maximum NA, here 1.49 1^ . 
For Fourier analysis, a 50/50 mirror was used instead of the dichroic beamsplitter. 

Numerical simulations were realized with a commercial Finite-Element-Method software. 
The full three-dimensional model was composed of an orthorhombic gold nanorod placed on 
a glass substrate with a dielectric constant £^=2.25. The gold dielectric constants were taken 
from l(30i . Perfectly matched layers were added to the model to avoid any boundary refiections. 
The model was solved for an incident plane wave coming from the substrate with a propa¬ 
gating direction orthogonal to the glass surface. To take into account the fast decrease of the 
electromagnetic field around the nanorod, the discretization was finer than 2 nm. 

Results and Discussion 

Figure [^a) shows the experimental TPF intensity map as a function of nanorod length obtained 
for a HGo^oi^^y) excitation beam. The horizontal axis corresponds to the spatial position of the 
nanorods relatively to the focused beam center. Because TPF intensity has a quadratic depen¬ 
dence on the near-field intensity an increase signal relates to a local enhancement of the electric 
field; TPF measurements are thus an efficient tool to identify FSP resonances ifTTll . The repre¬ 
sentation clearly shows an increase of the TPF intensity related to SP resonances for lengths 
comprised in the [95,130] nm and [330,400] nm ranges. Figure [^b) reports two enlarged con- 
focal maps for 96 and 780 nm long nanorods. Confocal maps result from the convolution of the 
spatial intensity distributions of the excitation beam and of the SP mode. These maps are far- 
field images and their resolution is then diffraction-limited to 270 nm. Thus, nanorods shorter 
than this latter length directly probe the excitation intensity profile. In view of the confocal map 
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Fig. 2. Bright Modes : (a) Distribution of TPL intensity versus nanorod length for a 
//Go,o (-^5 j) excitation beam. The horizontal axis corresponds to the nanorod position rel¬ 
atively to the excitation beam, (b) Zoom on the confocal images of the TPL from the 96 
and 780 nm long nanorods. The scale bar is of 500 nm. (c) The length dependence of the 
maximum TPL signal (red columns) is compared to the simulated intensity of the electric 
field integrated 5 nm around the nanorods (blue line) for a plane wave excitation, (d) Simu¬ 
lated spatial distribution of the Ex component 10 nm below the nanorod for a 100 nm and a 
380 nm long nanorods at an excitation wavelength of 800 nm. The black lines indicate the 
contour of the corresponding nanorods. 


of the 96 nm long nanorod, the excitation appears to have a gaussian profile as expected for 
the related excitation conditions. For longer nanorods, the TPL is only observed at the nanorod 
extremities and results from a local enhancement of the electric field at the structure boundaries 
as shown for the 780 nm long nanorod. 

In Fig.|^a), only the length dependence of the TPL mapping is reported and Fig.|^c) reports 
the length dependence of the maximum of TPL intensity. The first resonance occurs for length 
of 96 nm where the TPL signal is 25-fold higher than the one of off-resonance nanorods. The 
experimental data (red columns) are compared to simulations (blue solid line). Calculations 
of the near-field of orthorhombic structures produce strong field enhancements at the corners 
due to the abruptness of the structure. In order to reduce the impact of this field enhancement 
compared to the exaltation due to surface plasmon resonances, we report here the mean electric 
field intensity integrated in a volume of 5 nm around the nanorods instead of its maximum. 
Simulations are realized with a plane wave illumination linearly polarized along the nanorods. 
Calculated plasmonic resonant lengths are in good agreement with those experimentally de¬ 
tected by TPL measurements. Then, we plot in Figure [^d) the distribution of the simulated 
Ex component (oriented along the nano wire long axis) for the two resonant lengths: 100 and 
380 nm. These maps show that, at resonances. Ex is maximum at the center of the nanorod 
which ensures a good coupling with a focused l^ser beam once it is centered to the 

nanorod. Here, we precise that these near-field Ex maps cannot be compared to confocal TPL 
maps. To be convinced, one just has to note that simulated near-field Ex components oscillate on 
distances shorter than 130 nm whereas TPL maps have a diffraction-limited resolution of 270 
nm since they are far-field images related to E^. The resonance lengths (100 and 380 nm) corre¬ 
spond to the first (m=l) and third (m=3) order plasmonic modes since Ex presents respectively 




































one and three field maxima inside the nanorods (the m index is the ratio between the nanorod 
length and the effective wavelength as described below in the Fabry Perot description of reso¬ 
nances). However, as a change of sign in the distribution corresponds to a ;r-phase change, 
these odd modes present two ;r-phase changes at their extremities, as previously reported ED 
The Coupled Mode Theory proposed by A. Yariv 1i32\ says that optical modes are efficiently 
excited when the excitation intensity and phase spatial profiles match with those of the latter 
optical modes. Then, we believe that SP modes could be more efficiently excited with focused 
HG2^o (-^5 y) or HG4^o (-^5 y) beams for which the phase-matching would be optimized. Here, res¬ 
onances of longer nanorods with higher-order odd modes are not observed experimentally since 
the overlap with the HGo,o{x^y) is vanishing when the mode order (m) increases. 
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Fig. 3. Dark Modes : (a) Mapping of the nanorod TPL intensity versus the nanorod length 
for a //Gi o(v,y) excitation beam. The horizontal axis corresponds to the nanorod position 
relatively to the excitation beam, (b) Zoom on the confocal images of the TPL from the 
96 and 250 nm long nanorods. The scale bar is of 500 nm. (c) The length dependence 
of the maximal TPL intensity (red columns) is compared to the simulated intensity of the 
electric field integrated in a volume of 5 nm around the nanorods (blue line). The excitation 
is similar to a //Gi o(v,y) focused beam centered on the nanorod. (d) Simulated spatial 
distribution of the Ex component 10 nm below the nanorods for lengths of 250 nm and 
520 nm at an excitation wavelength of 800 nm. 


Next, we interrogate the TPL intensity for a HGi^o{x,y) excitation beam in order to excite 
even eigenmodes. Results are reported in Fig.which is organized in the same way than Fig.[^ 
Here, simulations are realized with an excitation similar to an HGi^o{^,y) focused beam. The 
TPL map of Fig. [^a) clearly shows a rise of the nonlinear response for a new set of length 
ranges: [240-300] nm and [480-580] nm. However, this figure is slightly more complex than 
Fig. [^a). Indeed, as previously mentioned, confocal maps result from the convolution of the 
spatial intensity distributions of the SP modes and of the excitation beam which is here modified 
by the spatial phase shaping. As above mentioned, shortest nanorods probe the spatial excita¬ 
tion intensity distribution. Figj^b) shows that the confocal TPL map of 96 nm long nanorod 
is mainly composed of two lobes as expected for a focused HGi^oi^^y) beam. Note that the 
asymmetric lobe intensities result from imperfections of the home-made A/4-step gold mirror. 
Enlarged confocal map of the resonant 250 nm long nanorod exhibits three TPL maxima. The 
most intense one is obtained when the laser beam is centered on the nanorod, i.e. when the 











































phase profiles of the excitation and of the plasmonic mode coincide. The two minor lateral 
spots are obtained when the nanorod extremities overlap each lobe of the focused HG\^{){x^y) 
beam. Even modes are thus efficiently excited only when the nanorod is centered with the fo¬ 
cused laser beam. For longer nanorods, confocal maps are becoming increasingly 

complex with four or even five lobes resulting from the convolution of confocal images simi¬ 
lar to Figj^a) with the ifGio(v,y) intensity distribution. The understanding of these complex 
maps is beyond the scope of this work. In the following, we consider the excitation of even 
modes occurring in the centered condition which is illustrated by the gray box of Fig. [^a). 
Figure, [^c) reports the maximal TPL intensity in that centered condition together with simu¬ 
lated near field intensities calculated with an excitation similar to a HG\^Q{x^y) profile. Here 
again, there is a good agreement between the calculated length dependence of the electric field 
intensity (blue line) and the experimental data (red columns). Figure, [^d) gives the simulated 
Ejc spatial distribution for both resonant lengths: 250 and 520 nm. These modes are the second 
(m=2) and fourth (m=4) order eigenmodes since 2 and 4 maxima of are respectively present 
inside the nanorod. The dark nature of these modes is confirmed by the presence of a ;r-phase 
change at the center of the nanorods which also ensures a good match with the HG\^Q{x^y) 
phase profile. Here, we unambiguously demonstrate the selective excitation of dark plasmonic 
modes by phase-shaping the excitation beam. One has to note that dark modes do not radiate 
by essence. Nevertheless, TPL which does not depend of the mode dipolar nature is still ex¬ 
alted at resonances (even if dark) since it is only sensitive to the local field intensity. TPL is 
then powerful to probe both bright and dark resonances of SPs. Finally, we anticipate that more 
complex phase-shaped beams (eventually produced by spatial light modulators) could be used 
to efficiently excite higher-order bright or dark modes. 

LSP resonances can be described by the commonly accepted Fabry Perot resonator descrip¬ 
tion of localized surface plasmon resonances EEiiiia. Here, the fundamental propagating 
mode of a gold nanowire (not shown here) with a cross-section of 55 nmx35 nm has a cal¬ 
culated effective index Ueff of 3.03. The effective wavelength Xeff is thus ^lueff=264 nm. 
The tail of the field at the nanorod ends decays with a characteristic distance close to 5=12 nm 
(resp 6, 4, 2.5 nm) for the m=l (resp. m=2, 3, 4) mode (see Fig|^d) and Fig|^d)). In the 
above-mentioned description, the resonant modes are occurring when L + 25 = 

Resonances are thus expected for lengths of L=108 nm, 252 nm, 388 nm, and 523 nm, in 
quite good agreement with the experimental four plasmonic resonances measured for lengths 
of L=100 nm, 250 nm, 380 nm and 520 nm. 

Now we turn to a far-field analysis of scattering properties of the modes by Fourier imaging. 
Bright and dark LSPs differ by their near-field distributions, and also by their far field scattering 
properties which are also dictated by the order of the modes ||34j[35l. In Fig.[^a-e) we report 
the experimental Fourier-plane images of the back-scattered light for the four resonant modes 
above defined and for an out-of-resonance 909 nm long nanorod. Detection is polarized along 
the Ex component (parallel to the nanorod long axis). In these images, the central black circle 
corresponds to the solid mask placed in the detection path way to reject the light distributed 
within the excitation numerical aperture. Figures (a), (b) and (c) concern the HGo^o{x^y) 
excitation of the m=l and m=3 bright LSP resonances and of the non resonant nanorod. These 
Fourier-plane images show that the bright modes principally radiate in a direction orthogonal to 
the nanorod long axis as expected from a dipolar emission. In particular, shortest nanorods have 
a dipolar-like pattern whereas the longest one exhibit fringes due to interferences originating 
from light scattered from different nanorod positions. It results from these interferences that the 
central fringe is narrower for longer nanorods and the directivity is then increased. 

Following the same approach, the images of Fig. [^(d) , (e) and (f) concern the HGi^oi^^y) 
excitation of the m=2 and m=4 dark LSP resonances and from the non resonant nanorod. It is 
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Fig. 4. (a-e) angular distributions (Fourier plane) of the back-scattered light polarized 
along the nanorod for m=l to m=A modes, and for an off-resonant 909 nm long structure. 
The detection is limited by the NA of the objective (1.49). The critical angle at NA=1 is also 
marked. The filtering mask in clearly observed at the center of the Fourier plane images, (g- 
1) angular distributions calculated with Dyadic Green function for above mentioned cases, 
(m-r) angular distributions calculated by the coherent superposition of the field emitted by 
two in-phase or TT-dephased dipoles separated by the corresponding nanorod lengths. 


worth to note that dark modes poorly scatter in far-field, and the experimental signal from these 
modes is at least ten times lower than bright modes. More precisely, for antenna shorter than 
500 nm, we only collect signal from resonant antenna. The Fourier images of Fig. (d) , (e) 
and (f) clearly differ from those of bright modes. In particular, dark modes present a minimum 
of emission exactly where bright modes principally emit. The orthogonality of the local field 
distribution between these modes appears here: in this plane orthogonal to the nanorod axis, 
bright modes constructively interfere whereas dark modes destructively interfere. Indeed, the 
HG\^Q{x^y) excitation generates out-of-phase plasmon in both nanorod sides. Furthermore, it 
is especially interesting to note that the directivity of the emission from 250 nm long resonant 
nanorod is tilted by 90 degrees from the directivity of brights modes. Curto et al. reported sim¬ 
ilar emission patterns resulting from a near-field excitation of single nanorods with a quantum 
dot placed in the vicinity of plasmonic elements ll^ . In the present work, the radiation patterns 
are obtained by scattering the laser field onto the antennas. Since the same modes are excited, 
the angular symmetry of the scattered radiation is the same as in HSU- However, the origin of 
the signal is fundamentally different. With our approach, we demonstrate that the control of 
the spatial phase-profile of the optical excitation is a powerful tool to engineer the directivity 
of an optical antenna. Furthermore, from our measurements it seems that in the Fourier images 
do not conserve information on the mode orders and are only governed by the nanorod diffrac¬ 
tion figures. More precisely, resonant nanorods have emission patterns similar to off-resonance 
nanorods of similar lengths (not shown) and only the intensity is modified. Here, the Fourier 
images are only sensitive to the parity of the excited mode. 

Experimental emission patterns are also confronted to th ose calculated using the Green’s 
dyad technique. To this aim, we first calculate the electric field inside the nanostructure as 
described in 1361 . For simplicity, we use a plane wave excitation to simulate the HGQy){x^y) 
excitation and a plane wave excitation but with a ;r-phase shift for negative x to compare to 
the HG\^{){x^y) excitation. Then the intensity scattered in the Fourier plane is calculated using 
the asymptotic expansion of the Green’s dyad (371 [38l. Calculated Fourier-plane images are 
reported in Fig.|^(g-1). The main features of experimental images: 1) the increase in directivity 










for longest nanorods, 2) the change in orientation of emission from bright to dark modes and 
3) the zero emission in the nanorod perpendicular plane for dark modes are well described 
by the calculations. This technique also describes the fine structure of the emission patterns 
of the excited longest nanorods for which two lobes appear on both sides of the 

main fringe. If this technique is in good agreement with experimental measurements and allows 
to anticipate the emission properties of plasmonic nano-antenna, it is highly time consuming 
especially for multi-element and large scale antenna like Yagi-Uda. 

Then we propose a simple model to describe the nanorod emission pattern which is only 
based on the emission of two dipoles. In this basic model only two dipoles are considered. De¬ 
pending on the dark or bright mode, they oscillate in phase or they are ;r-dephased as depicted 
at the left of figure |^m) and (p). Dipoles are placed 10 nm above a silica substrate (n=1.5). 
Figure^m-o) (resp. (p-r)) report the calculated Fourier-plane images for in-phase (resp. out-of- 
phase) dipoles separated by the corresponding nanorod lengths. Qualitatively, this quite simple 
model reproduces very well the main features of the experimental patterns especially for short¬ 
est antenna. The accurate comparison with experimental emission patterns of bright and dark 
modes of short antenna is quite surprising and gives us a better understanding on the origin 
of the angular emission patterns. Indeed, the physical meaning of this two-dipole based model 
corresponds to a Fourier-plane pattern essentially governed by radiations from the nanorod ex¬ 
tremities. It indicates that the radiative losses along the nanorods can be neglected for antenna 
shorter than 700 nm and only extremities emissions can be considered. Then, bright or dark 
modes only differ from the phase between extremities emissions. For longer antenna, the radi¬ 
ations all along the nanorod start to take an important role in the scattering patterns and have to 
be taken into account. The resulting patterns are then governed by diffraction figures as reported 
by Sersic et al 1391 . However, we anticipate that the simplicity of this two-dipoles model can 
help to estimate the emission patterns of complex multi-element structures interacting through 
bright and/or dark modes. Furthermore, it has the advantage to propose analytical solutions and 
is poorly time-consuming compared to full-3D method like Green or FEM methods. 

Conclusion 

To summarize, we used tightly focused Hermit-Gaussian modes to selectively excite bright and 
dark plasmonic eigenmodes sustained by individual gold nanorods. Non-linear photolumines¬ 
cence was used as probe of the local field maxima generated by the even or odd natures of the 
modes. In particular, we were able to unambiguously demonstrate the excitation of dark modes 
by introducing a ;r-shift in the excitation beam. We further showed that the emission pattern 
features the symmetry of the mode and can be qualitatively understood by considering two 
dipoles oscillating in phase or out of phase at the nanorod extremities. This work provides a 
simple avenue for selectively exciting bright or dark resonances that can be subsequently used 
to control the plasmonic response of single optical antennas. 
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